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ABSTRACT 

Aims. We derive the age and metallicity of circumnuclear star formation regions (CNSFRs) located in the spiral galaxies NGC 695 1 and 
NGC 1097, and investigate the cause of the very low equivalent widths of emission lines found for these regions. 

Methods. We used optical two-dimensional spectroscopic data obtained with Gemini GMOS-IFUs and a grid of photoionization models to 
derive the the metallicities and ages of CNSFRs. 

Results. We find star formation rates in the range 0.002-0. 14 M Q yr and oxygen abundance of 12+log(0/H)ai8.8 dex, similar to those of most 
metal-rich nebulae located in the inner region of galactic disks. 

Conclusions. We conclude that the very low emission-line equivalent widths observed in CNSFRs are caused by the "contamination" of the 
continuum by (1) contribution of the underlying bulge continuum combined with (2) contribution from previous episodes of star formation at 
the CNSFRs. 
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1. Introduction 



The pioneer works by Morgan (1958) and Sersic & Pastoriza 
(1967) have shown that many spiral galaxies host circumnu- 
clear star formation regions (CNSFRs) and several works have 
studied these regions. In one of the most recent works, Knapen 
(2005 ) uses Ha images of a sample of spiral galaxies to study 
CNSFRs, as well as nuclear star formation regions. They found 
that CNSFRs are present in about 20% of spiral galaxies and al- 
most always occur within a barred host. In this paper we adopt 
the definition of Knapen (120051 ) to distinguish CNSFRs from 
nuclear star formation regions: CNSFRs are the ones observed 
within circular areas with diameters of a few kpc around the 
galactic nucleus, while nuclear star formation regions are ob- 
served in more internal areas a few hundred parsec across. In 
particular, the formation of CNSFRs in barred galaxies seems 
to be due to the radial inflow of gas along bars to the galac- 
tic center (Roberts et al. 11 9791 Friedli & Benz ll995l Sakamoto 
et al. [19991 Sheth et al. 120051 Jogge et al 120051 . This gas is 
accumulated in a ring structure between two inner Lindblad 
resonances (Knapen et a U1995t . 

Gas inflow is also necessary to feed the nuclear massive 
black hole in active galactic nuclei (AGNs) and surround- 
ing nuclear star formation regions, probably yielding the so- 
called AGN-starburst connection (Norman & Scoville 1988; 
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Terlevich et al. 1 19901 Heckman et al. 1997; Gonzalez Delgado 
et al. [19981 Ferrarese & Merritt 120001 Cid Fernandes et al. 
I200T1 Storchi-Bergmann l200Tl Heckman l2004b . Recent obser- 
vational studies have shown evidence of gas inflow: Storchi- 
Bergmann et al. ( 120071 1 and Fathi et al. (120061 ) report the dis- 
covery of gas streaming motions along spiral arms towards 
the LINER nuclei of the galaxies NGC 6951 and NGC 1097, 
while Fathi et al. (2007) found gas inflow to the nuclear star 
formation region of the galaxy M 83. Moreover, several works 
have found star formation around active nuclei (e.g. Riffel et al. 
120071 Kauffmann et al. 120031 

In CNSFRs, the line ratio [Om]A5001/H/3 (hereafter 
[Om]/HjS) is usually weak or not observable and the emis- 
sion lines have lower equivalent widths than those of disk H n 
regions (e.g. Diaz et al. 120071 Bresolin et al. 1999; Bresolin 
& Kennicutt 119971 Kennicutt et al. 119891) . The low value of 
[Om]/Hj6 has been attributed to over-solar abundances (e.g. 
Boer & Schulz [T9931 l. At their location in the innermost parts 
of galactic disks, these objects are expected to be the most 
metal-rich star-forming regions. In fact, Diaz et al. (2007) ob- 
tained long-slit observations for 12 CNSFRs and, using a semi- 
empirical method to derivate of oxygen abundances, found that 
some of them have 12+log(0/H)=8.85, a value consistent with 
the maximum oxygen abundance value derived for central parts 
of spiral galaxies (Pilyugin et al. l2007b . 

Regarding the equivalent width of emission lines, 
Kennicutt et al. ( 1989 ) found that the equivalent widths of Ha 
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Fig. 1. Sample spectra: star-forming regions F in NGC6951 
(top) and J in NGC 1097 (bottom). 

EW(Ha) in CNSFRs are about 7 times lower than the ones 
in disk Hn regions. They proposed several scenarios to ex- 
plain this behavior, such as (i) deficiency of high-mass stars 
in the initial mass function (IMF), (ii) a long timescale for 
star formation, and (iii) high dust abundance. However, nei- 
ther of these mechanisms turned out to be the dominant effect. 
Unfortunately, spectroscopic data and abundance estimates are 
only available for a small number of CNSFRs. 

With the goal of increasing the number of CNSFRs with 
determined gas abundances and star formation rates, we have 
combined two-dimensional integral field unit (IFU) data of 
CNSFRs located in the galaxies NGC 1097 and NGC 695 1 with 
a grid of photoionization models to derive their nebular gas 
properties. Using high-quality two-dimensional data, we could 
also investigate the cause of lower equivalent widths of emis- 
sion lines observed in CNSFRs with respect to those of inner 
disk HII regions. In Section[2]we give details about the obser- 
vational data we used, in Section[3]we present a description of 
the methodology used to derive the physical parameters, while 
in Section [4] we present our results. The discussion of the out- 
come and our final conclusions are given in Sections and [6] 
respectively. 

2. Observational data 

The observational spectroscopic data on NGC 1097 and 
NGC 6951 were drawn from the works of Fathi et al. (2006) 
and Storchi-Bergmann et al. (2007). Basically these data were 
obtained with IFUs of the Gemini Multi-Object Spectrographs 
(GMOS-IFU) on the Gemini South (for NGC 1097) and North 
(for NGC 6951) telescopes using the R400-G5325 grating and 
r-G0326 filter, covering the spectral range 5600-7000 A with 
a resolution of R * 2300. Three IFU fields were observed in 
each galaxy, covering an angular field of 15"x 7". The reduc- 
tion procedure resulted in 50 x 70 spectra for each IFU field, 
each corresponding to an angular coverage of 0.1"x 0.1". We 
adopted distances for NGC 695 1 and NGC 1097 of 24 Mpc and 
17 Mpc (Tully 1988), respectively, such that l"corresponds to 



96 pc for the former and 68 pc for the latter galaxy. The reader 
is referred to the papers above for a full description of the data. 

Figure Q] shows spectra of CNSFRs identified by the let- 
ters F and J in Figs.[2]and[3]located in the circumnuclear rings 
of NGC 6951 and NGC 1097, respectively. Using the observed 
spectra we built maps for flux distributions in the Ho- emis- 
sion line and for the line ratios [Nn]/t6583/Ha, ([Sn]/16717 + 
A6731)/Ha, as well as for the EW(Ha). These maps are shown 
in Figs. [2] and [3] Throughout the paper we identify the line ra- 
tios above as [Nn]/Ha and [Sn]/Ha. For NGC 6951 the ob- 
served field completely covers the circumnuclear ring. The data 
of the ring were separated from those of the inner region by fit- 
ting two ellipses to the Ha image of NGC 695 1 : one to the inner 
border of the ring and the other to the outer border (see Storchi- 
Bergmann et al. 120071 1. In the NGC 6951 Ha map, we clearly 
see several regions of star formation. We have identified nine of 
them as indicated in Fig. [2] Since for NGC 1097 the observa- 
tion field covered only part of the ring, we could only identify 
two regions in this galaxy (see Fig. 0). 

For these measurements, we assumed that each star forma- 
tion region has a circular symmetry, whose center corresponds 
to the peak in the Ha intensity and the radius defined as the one 
where the flux reaches 50 % of the peak value. Nevertheless, we 
can see from Fig.|2]that the CNSFRs have irregular morpholo- 
gies and the assumption above can result in some uncertainties 
in the Ha flux and in line ratio measurements. To estimate how 
much flux we may be missing, for NGC 6951, we compared 
the sum of the Ha flux of all regions (obtained from photom- 
etry within circular apertures) with the flux integrated over the 
whole ring. We found that the diffuse gas contribution (from 
regions not covered by the circles) is about 12 %. A similar re- 
sult was found for individual H n regions by Oey & Kennicutt 
(119971) . In relation to line ratios, we measured the median val- 
ues considering different values for the radius of each region 
and differences no greater than 10% were found. We thus es- 
timate that uncertainties in both fluxes and line ratios due to 
uncertainties in the geometry of the CNSFRs are about 10%. 
In Table [T] we show the adopted area of each region, the lumi- 
nosity of Ha, and the median values of EW(Ha), [N n]/Ha, and 
[S n]/Ha. 

3. Photoionization models 

To derive the metallicity, ionization parameter, and the age 
of the CNSFRs, we employed the photoionization code 
Cloudy/95.03 (Ferland 2002) to build a grid of models cov- 
ering a large space of nebular parameters. The models were 
built as in Dors & Copetti (120061 ). with metallicities Z= 1 . 1 , 1 .0, 
0.7, 0.6 Zq, and ionization parameter logarithm (log U) of -2.0, 
-2.5, and -3.0 dex. In each model a stellar cluster was assumed 
as the ionizing source with the stellar energy distributions ob- 
tained from the synthesis code STARBURST99 (Leitherer et 
al.[l999). We built models with stellar clusters formed by in- 
stantaneous and continuous bursts with stellar upper mass lim- 
its of M up = 20, 30 and 100 M , Z= 1.0 and 0.4 Z Q , and ages 
ranging from 0.01 to 10 Myr with a step of 0.5 Myr. Clusters 
older than about 10 Myr and formed instantaneously have no 
more O stars (Leitherer & Heckman 1995), so were not consid- 
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Fig. 2. Star-forming ring in NGC6951. From left to right: map of the Ha integrated flux (ergs cm 2 per pixel); line ratio maps 
[Nul/Ha [S n]/Ha and EW(Ha). The circles and the letters identify each region as defined in Section|2] 



ered in our models. The metallicity of the nebula was matched 
to the closest available metallicity of the stellar cluster. The so- 
lar composition refers to Asplund et al. (2005 ) and corresponds 
tolog(0/H) = -3.30. 

4. Results 

Since the observational data for NGC 695 1 and NGC 1097 cov- 
ers a restricted wavelength range, few emission line ratios can 
be used in our analysis. Storchi-Bergmann et al. ([1994) and 
Raimann et al. (2000 ) showed that [N ii]/Ho' can be used to ob- 
tain the oxygen abundance of star-forming galaxies. However, 
this line ratio also dependents on the ionization parameter U 
(Kewley & Dopita 2002 ), thus the combination with other line 
ratios is desirable in order to decrease the uncertainties in the 
O/H estimates. Thus we used our photoionization model grid 
to build an [Nn]/Hor vs. [S n]/ Ha diagram to derive Z and U 
of the CNSFRs located in NGC 6951 and NGC 1097. Recently, 
Viironen et al. (2007) have shown that diagnostic diagrams us- 
ing [N n]/Ho' and [S n]/ Ha are a powerful tool for estimating 
Z and U of star-forming regions. Moreover, these line ratios 
are not dependent on the M up and age assumed in the models 
(see Fig. [5]). We point out that a cautionary note has been put 
forth by Mazzuca et al. (2006) regarding the use of diagnostic 
diagrams including the line ratio [N ii]/Hq': they can yield de- 
generate values for the metallicity, as star-forming regions with 
low Z and U have [N n]/Ha values similar to regions with high 
Z and U. However, CNSFRs are known to have high abun- 
dances (near solar) (see Sect. 15.11 ). ruling out the possibility 
of low abundances. This is confirmed by the work of Storchi- 
Bergmann et al. {1996 ), which shows an increasing gas abun- 
dance towards the nucleus, including a few regions in circum- 
nuclear rings. 

The results of our models compared to observational data 
for NGC 6951 and NGC 1097 are shown in Fig. |U in which 
we considered an instantaneous burst with age of 2.5 Myr and 
M m = 100 M Q . We can see that the observational data are 



well represented by models having metallicities Z * Z Q and 
log U * -2.3, with uncertainties of about 0.03 and 0.1 dex, due 
to the uncertainties in the line ratios discussed above. In TableQ] 
we show the Z and log U values for each region, obtained by 
linear interpolation from the model grid. In this table we also 
present the values of electron density N e computed from the 
observed line ratio [S ii]/16717//i6731 and adopting T e = 10 4 
K (Osterbrock 1989). The N e values cover the range 50-600 
cm" 3 , which is typical of Hn regions (Copetti et al. 2000). 

To obtain the ages of the CNSFRs, we used the observed 
values of equivalent widths and [O m]/H/3, which decrease al- 
most monotonically as a function of time (Magris et al. 2003 
Copetti et al. 1986 ). In Fig.|5]we plot the values predicted by the 
models for EW(Bry), EW(H/3), EW(Ha), as well as [Nn]/Ha, 
[S n]/Ha and [O m]/H/3 versus age. We considered instanta- 
neous and continuous star formation bursts with M up = 20, 30 
and 100 M , log U - -2.5 and Z = Z . The values of the line 
ratios [N n]/Ha and [S n]/Ha were found to be approximately 
constant for the range of ages considered. However, [O m]/H/3 
shows a chaotic behavior between 2 and 6 Myr, mainly for 
M up =100 M and instantaneous burst. This happens due to the 
increase in the number of ionizing photons by the presence of 
Wolf Rayet stars in this phase (Leitherer & Heckman 1995). 
This behavior was also noted by Dopita et al. (2006). As no 
feature of Wolf Rayet stars was observed in our data and be- 
cause [O m]/HyS is very dependent on many nebular parame- 
ters, it will not be considered as an age tracer in this paper. The 
use of the equivalent width of emission lines seems to be more 
reliable for deriving ages of star-forming regions, since they 
basically depend on M up and on the slope of the IMF (Copetti 
etal. [TM6i ). 

Using the EW(Ha) values shown in Table Q] the CNSFR 
ages predicted by instantaneous and continuous models are 
around 7 Myr and over 10 Myr, respectively, independent of 
the M up assumed. Table Q] also presents the values for the star 
formation rate (SFR) computed using the luminosity of the re- 
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Fig. 3. Left: map of the Ha-integrated flux (ergs crrT 2 per pixel) showing the field covered by the IFU in NGC 1097, which 
includes the nucleus and part of the ring (at the bottom). Right: line ratio maps [Nn]/Ha [S n]/Ha and EW(Ha) for the partial 
ring. The circles and the letters identify each region defined in Section[2] 
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Fig. 4. Diagnostic diagram [Nn]/Ha vs. [Sul/Ha showing the 
grid of photoionization models for an ionizing cluster with age 
2.5 Myr and M up = 100 M B . Solid lines connect curves of iso- 
Z, while dotted lines connect curves of iso-f/. The values of 
log U and Z are indicated. Squares are the median values mea- 
sured for the CNSFRs in NGC 695 1 and triangles are those for 
NGC 1097. 
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Fig. 5. Evolution of the EW(Bry), EW(H/3), EW(Ha), and rel- 
evant line ratios as predicted by our photoionization models 
considering different M up values and star formation regimes as 
indicated, for Z = Z Q and log U = -2.5. 



gions obtained in this paper and the relation SFR(M Q year ')= 
7.9 x L(Ua) (ergs s~') given by Kennicutt dl9981 >. 



5. Discussion 

5.1. Metallicities 

Acording to our models, the CNSFRs of NGC 6951 and 
NGC 1097 have solar or slightly higher solar metallicities. 
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Table 1. Physical properties of the CNSFRs in NGC6951 and NGC 1097 
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Considering an uncertainty on the order of 0. 1 dex due to the 
method used to estimate abundances (see Evans 1986, Dors & 
Copetti 2005 , 2006 ), our abundance values are comparable to 
the maximum oxygen abundance derived for central parts of 
spiral galaxies of 12+log(0/H) * 8.87 (Pilyugin et al. l2"0"07j l. 
Similarly, Mazzuca et al. (2006) present an emission-line diag- 
nostic analysis of the nuclear starburst ring of NGC 7742 and 
find roughly solar abundances. 

Recently, Sarzi et al. ( 120071 1 presented a study of CNSFRs 
located in eight galaxies which include the regions B, E and F 
of our sample. They also used empirical methods to derive the 
metallicities of these regions and find abundances similar to 
the ones derived in this paper. Diaz et al. (2007) use a semi- 
empirical abundance calibration to find that the CNSFRs in 
NGC 2903, NGC 3351, and NGC 3504 have metallicities com- 
parable to the ones derived in this paper. Until now, only a 
case of very low abundance on the order of 0.2-0.4 Z Q in the 
CNSFRs located in NGC 3310 was reported by Pastoriza et al. 
(1993 ). This low abundance can be understood as due to the 
fact that this galaxy has an interaction history (Elmegreen et al. 
2002 ), through which neutral gas may have been accreted to 
the ring. 

Due to the fact that CNSFRs have very low excitation, with 
[Om]/HjS < 1 (e.g. Hagale et al. [20071 Perez et al. [2000b 
and to their central location in the galactic disk, it is expected 
that these objects have somewhat higher abundance than is de- 
rived for inner disk Hn regions. To check this, we have plot- 
ted the [Nn]/Ha in Fig. [6] (lower panel) against oxygen abun- 
dance for some disk Hn regions located in 13 spiral galax- 
ies. This sample includes the most metal-rich nebulae in which 
temperature sensitive emission lines have been detected. The 
data were taken from Bresolin (2007), Bresolin et al. (|2005 
120041 ), Kennicutt et al. ([20031 , Castellanos et al. d2002b . and 
Skillman et al. ( 119961) . With the exception of the oxygen abun- 
dances from Skillman et al. ( 119961 ). which were recalculated 
by us using the P-method (Pilyugin 2001 ), all O/H values were 
obtained using direct derivation of the electron temperature. 
Also included in Fig. [6] are the data of our sample. We can 
see that the CNSFRs show the highest oxygen abundances, to- 
gether with some disk nebulae, a result also obtained by Diaz 
et al. (2007). Since abundance determinations via photoioniza- 
tion models are overestimated by a factor of 0.1-0.4 dex when 



compared to direct determinations (Dors & Copetti 2005 ), our 
abundance estimates should be interpreted as upper limits. If 
CNSFRs are included in the computation of abundance gra- 
dients in spirals, we would expect therefore a plateau at very 
small galactocentric distances. 

5.2. Star formation history 

The Ha luminosity values found for the CNSFRs in this paper 
are in the range 2 x 10 38 - 1.2 x 10 40 erg s _1 , compatible with 
values of regions in circumnuclear rings of other spiral galaxies 
(e.g. Diaz et al. 2007). These luminosity values yield star for- 
mation rates of 0.002-0.1 M yr _1 , characterizing a moderate 
star-forming regime. Averaging the SFR values in Table [TJ we 
find 0.04 M yr _1 , a similar value to the one found by Ho et 
al. d 19971 ) for a sample of nuclear star-forming regions and by 
Shi et al. (2006), who studied a sample of 385 circumnuclear 
star-forming regions in galaxies with different Hubble types. 

We cannot conclude from our data alone whether the star 
formation has been continuous or instantaneous in the rings, 
although the former scenario has not been favored by previous 
studies. For example, Sarzi et al. (2007) use absorption-line in- 
dex diagrams to show that the hypothesis of a constant star for- 
mation in their sample can be ruled out. Allard et al. (2006 ) find 
that the star formation in the ring of M 100 have occurred in a 
series of short bursts for the past 500 Myr or so. These results 
can be understood within the following scenario. Once mas- 
sive stars are formed, they cause turbulence and heating in the 
molecular cloud, inhibiting further star formation (Hartmann et 
al. 1200 II Blitz & Shu |198til ). Since no systematic age sequence 
is observed along the NGC 695 1 ring, the star formation mode 
seems to be as predicted by the popcorn model (Boker et al. 
2008 ), in which the entire ring begins to form stars at about the 
same time or at random times. A similar result was also derived 
by Mazzuca et al. (2008 ), who used photometric Ha imaging 
of 22 nuclear rings to find that only a very small fraction of 
them show age sequences along the ring. 

5.3. Age estimates 

Our modeling gives ages of about 7.0 Myr for an instanta- 
neous burst and more than 10 Myr for a continuous burst, with 
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the CNSFRs in NGC 1097 slightly older than in NGC6951. 
To compare our age estimates with those of other CNSFRs, 
we collected equivalent widths of more CNSFRs from the 
literature and estimated their ages using Fig. [5] The data 
were extracted from Diaz et al. (l20"07j l for NGC 2903 and 
NG C 3505 , Hagele et al. d2007l for NGC 3351, Allard et 
al. d2006b for M 100, Reunanen et al. (|2000] | for NGC 7771, 
and Wakamatsu & Nishida (TT980b for NGC 43 14, which give 
EW(H J 8), EW(Ha-), andEW(Bry) in the range 1-17, 10-40, and 
5-60 A, respectively. The predicted ages for these EW values 
are of 5-8 Myr and more than 10 Myr for instantaneous and 
continuous bursts, respectively, thus similar to the ones we de- 
rived for the CNSFRs in the present paper. Somewhat higher 
EW values have been obtained by Mazzuca et al. d20081 > for 
HII regions in nuclear rings of more than 20 galaxies. On the 
other hand, studies also utilizing photoionization models give 
ages of 1.5-3 Myr for Hn regions located in galactic disks and 
for H II galaxies (Dors & Copetti 120061 Stasinska & Izotov 
120031 Bresolin et al. [19991 Copetti et al. [19851 . We have also 
compared the ages of our sample with those of disk nebulae 
of similar metallicity. In Fig. [6] we show log(EW(Hy8)) and 
12+log(0/H) plotted against log([Nn]/Ha) for these disk Hn 
regions, together with our data for the CNSFRs. We converted 
the EW(Ha) values shown in Table Q] in EW(HyS) using the re- 
lation EW(Hy8)« 0.15 EW(Ha) obtained from our models. We 
can see that the EW(HyS) of the CNSFRs are lower by a factor 
of about 30 than those of their counterpart disk Hn regions. 
Assuming a EW(H/?)= 100 A for the innermost disk regions 
and using Fig. [5] we find that these regions are about 3 Myr 
younger than the CNSFRs. 

5.4. The low EW values 

The age differences found above may not be real, since the very 
low equivalent widths in CNSFRs stem from other effects than 
the age. Kennicutt et al. ([1989) pointed out possible mecha- 
nisms such as (i) deficiency of high-mass stars in the initial 
mass function, (ii) very high dust abundance, (iii) contamina- 
tion of the continuum by contribution from the bulge or other 
underlying stellar populations, and (iv) continuous star forma- 
tion. These hypotheses are discussed below. 

Regarding the hypothesis of the deficiency of high-mass 
stars, since M up of an ionizing cluster is dependent on the 
metallicity (Larson & Starrfield [192U KahnESll Shields & 
Tinsle yQl76l Stasin ska[l9 80l Vilchez & Pa gel[[988l Campbel l 
[19881 Bresolin et al. [19991 Dors & Copetti [20031 120531 l2006l l. 
CNSFRs and innermost disk nebulae should have about same 
M up . It is unlikely that the discrepancy in EWs is only due to 
variation in M up . Even assuming that the ionizing cluster of 
disk nebulae have M up = 100 M and the ones of CNSFRs 
have M up = 20 M , the age difference remains. Also a varia- 
tion in the IMF can be discarded, since stellar clusters appear to 
form with a universal IMF slope (Rroupa 2007, 2002, Chabrier 
120031 . 

The hypothesis of very high dust abundance in CNSFRs is 
not favored because they have similar Balmer decrements to 
those observed for disk H n regions, suggesting that this is not 



the main source of the low EWs in CNSFRs (Kennicutt et al. 
[19891 . 

Concerning the contamination of the continuum, several 
works have found the presence of underlying stellar popula- 
tions from previous generations of stars in CNSFRs, which can 
decrease the equivalent widths and yield no real age estimates 
from Fig. [5] For example, Buta et al. (2000) found a spread in 
ages of about 5 to 200 Myr for clusters in the circumnuclear 
ring of NGC 1326, with the largest number of clusters younger 
than 20 Myr. In the case of NGC 6951 and with measurements 
of the equivalent width of calcium triplet at a;8500 A, Perez 
et al. (2000) found evidence of a population of red supergiant 
stars with ages of 10 to 20 Myr in the ring. 

An underlying stellar population can contribute in two ways 
to the low EWs by causing a Balmer absorption and yielding 
an extra contribution to the continuum. To estimate the effect 
of an underlying Balmer absortion, we constructed models of 
stellar populations (using the Padova 1994 tracks) with ages in 
the range 10-1000 Myr. We normalized the continuum of each 
model to the observed one and subtracted them from the ob- 
served spectra of the CNSFRs of NGC 6951. For all cases, the 
increase in the EW values due to subtracting of an underlying 
absorption was at most 10%, much less than needed to com- 
pensate for the low observed EW in CNSFRs. 

Determining of the continuum contribution from underly- 
ing stellar populations is not an easy task, since two types of 
stellar populations can be contributing: the stellar population 
from the bulge and stars formed in previous episodes of star 
formation in the ring. To evaluate the contribution of the un- 
derlying bulge, we performed aperture photometry on the con- 
tinuum images (not shown) measuring both the CNSFR con- 
tinuum and that of the bulge in the surrounding regions. We 
concluded that the bulge contributes with « 75 % of the total 
flux in the case of NGC 695 1 and with 50 % in NGC 1097. The 
corresponding corrections are an increase by a factor of 2 and 4 
in EW(Ha) for NGC 1097 and NGC 6951, respectively. These 
factors increase the EW(Hcc) in NGC 695 1 from average values 
of*30A to ^120 A and from ^ 10 A to ^20 A in NGC 1097. 
The corresponding changes in age are from «7 Myr to *6 Myr 
in NGC 695 1 and no change for NGC 1097 («7 Myr). Thus the 
contribution from the bulge to the continuum also does not ex- 
plain the low EW's in CNSFRs. Stauffer ( fTWTb also studied 
the bulge contamination using aperture photometry of galaxy 
nuclei to measure the background contribution to the contin- 
uum flux, and Kennicutt et al. (119891 ) used those results to find 
a correction factor similar to the one found by us. 

Several generations of stars may coexist in CNSFRs be- 
cause the tidal effects are not strong enough to disrupt the clus- 
ter at the radius of the ring, allowing several generations of stars 
to coexist (Buta et al.[2000). Thus the continuum may have a 
contribution from previous bursts of star formation in the rings. 
To investigate the effect of this continuum, we built photoion- 
ization models to represent a scenario where three bursts of star 
formation have occurred, with ages of 0.01, 10, and 20 Myr for 
two values of M up : 100 and 30 M . This is illustrated in Fig. [7] 
which shows the evolution of EW(HyS) as a function of time. 
It can be seen that the EW values decrease considerably when 
compared with those of a single burst (Fig. [5]). If we consider a 
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Fig. 6. Oxygen abundance vs. [Nn]/Ha (lower panel) and 
EW(HyS) vs. [Nn]/Har (upper panel). The squares represent data 
collected from the literature, while the circles represent the data 
of the CNSFRs analyzed in this paper. 

typical EW(H/J) for CNSFRs of 10 A (not corrected for the un- 
derlying bulge contribution) and compute the age from Fig. [7] 
we find ages of 6-7 Myr, similar to the values derived above. 

We conclude that none of the mechanisms discussed above 
alone yield the low EW values (and resulting larger age) 
observed in CNSFRs when compared to disk Hn regions. 
However, it could be that a combination of more than one 
mechanism can explain the lower EWs in CNSFRs. For exam- 
ple, we can consider the effects of both the dilution by the bulge 
and by previous bursts of star formation. A typical EW(Hy8)= 
10 A in NGC 695 1 will increase to 40 A if we correct for con- 
tamination by the bulge. Then, we can consider the effect of 
previous generation of stars in the ring using Fig. [7] instead of 
Fig.[5]to obtain the age. We obtain an age of 3-4 Myr, a value in 
consonance with ages of the innermost disk regions. In this way 
we eliminate the apparent difference in age between CNSFRs 
and inner-disk H n regions. 

6. Conclusions 

We have used optical IFU spectroscopic data of CNSFRs in 
NGC 6951 and NGC 1097, which compared to a grid of pho- 
toionization models, in order to obtain physical parameters of 
these regions. We find that the star formation rates are in the 
range 0.002-0. 14 M yr _1 , thus covering a range from very low 
to moderate star-forming regime. We also find that the CNSFRs 
have oxygen abundances of 12+log(0/H)»;8.8, similar to those 
of the most metal-rich nebulae located in inner parts the disks 
of spiral galaxies. 

We investigated the cause of the very low equivalent width 
of emission lines found in CNSFRs, which suggests an older 
age for these objects than for disk H n regions. Among the sev- 
eral scenarios invoked to explain this fact, we have concluded 
that the contamination of the continua of CNSFRs by under- 
lying contributions from both old bulge stars and stars formed 
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Fig. 7. Evolution of the EW(Hy6) as predicted by photoioniza- 
tion models considering three episodes of star formation as in- 
dicated. The ages correspond to the youngest burst. 

in the ring in previous episodes of star formation (10-20 Myr) 
yield the observed low equivalent widths. Correcting for these 
contributions, there are not any more significant differences in 
ages between the CNSFRs and inner disk Hn regions. 
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